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TJw Motion of a Particle Attracted Towards a Fixed 

Center hy a Force Varying Inversely as the 

Fifth Power of the Distance. 

By William Duncan MacMillan. 



Introduction. 



In his " Traits des Fonctions Mliptiqiies," Legendre discussed briefly a 
number of central forces, one of which was the very general one 

In an article, published in 1853, Stader* considered a number of others in 
rather more detail, among them the force 

/=^ 

which has also been treated more recently by Miss Van Benschoten.f A short 
discussion of this law is given in many of the standard text books on Mechanics 
or Dynamics. The elliptic functions used in these discussions have been those 
of Jacobi and the treatment has been almost exclusively confined to orbits 
having apses. 

In the present paper we shall make no restrictions upon the initial con- 
ditions, except that they shall be real. The different types of orbits are exhibited 
as a one-parameter set of curves. It is shown that for every orbit having a real 
apse there exists another orbit such that the radius vector of one is proportional 
to the reciprocal of the other, for all values of the true anomaly, and that on 

»"De orbitis et motitius ptinotl culusdam corporei circa centrum attractionum aliis, quam Newtonia, 
attractlonls legibus sollicitati," Journal fiir Mathemaiik, Vol. XLVI, p. 862. 
t Master's Thesis in the Library of the University of Chicago. 
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every orbit not having a real apse there exists a real point such that if the true 
anomaly be measured from this point the radius vector corresponding to a positive 
value of the anomaly is proportional to the reciprocal of the radius vector for the 
negative value of the same anomaly. There is developed also a relation between 
the times which is analogous to Kepler's Harmonic Law. These relations are 
established in a very direct and elegant manner by the use of the elliptic functions 
of Weierstrass, which we shall, consequently, adopt. 

The paper is divided into two parts : In part I the problem is studied from 
the standpoint of the theory of functions without regard to reality questions ; 
in part II the cases of real orbits are discussed. 

Part I. 
§ 1. Differential Eq\mtion of the Orbits. 

For the force, /= -§, the diflferential equations in polar coordinates are 

where ^ is a constant depending upon the units chosen and is positive or negative 
according as the force is attractive or repulsive. Equation (2) furnishes at once 
the integral of areas 

^f = .. (3) 

We also find readily the vis viva integral 

From the relations 

dr dr dQ h dr 

df ~ dddi ~? dd 
we obtain 



c 



deJ 2h' ^ h^ ' ^^) 
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For arbitrary initial conditions 



d'i 



6=0, r = r„ -^° = a, (6) 



this equation becomes 

3 



\ dd/ ~ 2h'7i \ro) + [^ + '"' 2h'7i\ \ro) 



or, putting 

— 2 "> 



2AVg 

1 + . 
we get the equation 



' + "'-^1 = ^1^' 



0) 



If now we make the transformation 



(8) 



_ 1 _ 1 

we find 

/d^X 

l-l/=*''-(0+»(i)' « 

Equation (9) differs from equation (8) only in that b and /? are interchanged. 
Consequently any solution for — is also a solution for its reciprocal when b 
and jS are interchanged. 

§ 2. General Solution of Differential Equation. 

r 
Since only even powers of -7- occur in equation (8), let us put 



'0 

Prom this substitution there results 



{^)'=^^i'-Tr-T'^^-''0]['-j+j-^^^< (") 
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and this equation is reduced to the normal form of Weierstrass by putting 

ih 
whence ^ 

(^) = 4[« + i] \s-h-\^\-h^-] \_s-k + i^l-h^-]=R{s). (12) 
The roots of the equation R{s) = are 



e. = -i (13) 

e, = i_ JVl — 6/3; 
and 

eK + e^ + e, = 0. 

The solution of equation (12) is 

s = p{e + c"). 
Consequently 

The reciprocal relation gives us also the solution 

r^Y- ^^ ri5) 

UoV" p{e + c')-e, • ^^''^ 

In comparing these two solutions it is to be observed that e^,, e^ and e„ are 
unaltered by the interchange of h and /?, so that the p functions in the two 
cases have the same periods, these periods being functions of e^, e^ and e„. It is 
more convenient to take the reciprocals of (14) and (15), which are 

Taking the product of these two expressions we find 

ih^=[p{e + c")-e,-] lp{e + c')-e,-], 

and since this is true for all values of 6 it defines the relation between c' and c". 
If we take 6 =z — c" the first factor of the right member becomes infinite. 
Therefore the second factor must be zero, and consequently 

C' — C" = 6)^, 

where o^ is one of the half periods and p{a^) = e^. 
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If now we take 

c"=c— io^, 
our solution becomes 



and by virtue of the formula 

we have finally 



(17) 



(18) 



§3. Apses, Zero-Points and Infinity Points. 
Let us put 

0X = — C — i 6)^ + 6)x , 

^f. = — c—i 0)^ + 0)^, (19) 

e^ — —c — ^a^ + Oy, 

and denote the corresponding functional values by B^, J?^, By. We find then 

B, IT . V"P 

:^=0 =0, (20) 

To 

By IT . Vjb 



V- 



The derivative -=3- vanishes at the points 6 = 0x and 6 = 0^. These points 

are apses, while the point = B^ is a zero-point of the function. We may show 
this as follows : Put 
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Then (17) becomes . . 2 „ 

(y =|[«'(<^+"x)-a (21) 

and by means of the addition formula * 

and the formula 

p(r)-e,= [^^(r)]', 

it is readily reduced to 



r 



_ I 2 (e, - 6j f p (r) - e, 1* _ I 2 {e, - ej g, . . , . 



The (T quotient -^ (t) is an even function of -r; and therefore r(T) = r( — t), 
^\ 

and the point 6,^ is an apse. In a manner entirely similar it is shown that 0, is 
also an apse. For the point 6^, however, we find, by putting 6 = 6^ + t, 

(0' = |-b(^ + ".)-eJ; (23) 

and by the same formulsB as above this reduces to 

which is an odd function of t, vanishing for t" = 0. 
We have thus 

f(0, + r)= ^{O^-r), 

'0 '0 

L(e,+'r) = -^{d,-r), (25) 

'0 '0 

f(e. + '^)= f(e.-^). 

'0 ^0 

While the equations (17) have the periods 2a^, 2o^, 2o^, equations (18) 
have the periods 4ox, 4g)^ and 46),. f 
From the formula 

^{6 ± 2«.) = -^{6), ^(0 ± 2oJ =^(0), 

*See Schwarz, "i''orm«!« unci Lehrsatze eum Oebrauche der elUpiischen Fimctionen" §19, (5); also §18, (2). 
fSchwarz, §28, (3). 
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it is readily verified that the points 

0, + 26),, 0. + 2g),, 
0, + 26)^, 0. + 20)^, 

0;,+ 26),, 0, + 20), 

are all apses, while the points 

e^ + 26),, 
i + 2«,, 

0^ + 26), 

are all zero-points. In the complete period, therefore, there are 8 apses and 
4 zero-points. The points for which 6 equals 

— c — \a^, 

— c — la^ + 2(0;,, 

— G — \a^ + 2a>^, 

— c — \a^ + 26), 
are four infinity points. 

§ 4. Middle Points. 

Keturning to equations (18) let us take = — c and denote this value by 4'q 
and the corresponding value of r by B. If with these values we divide the first 
equation by the second, we find 

If we put = 4-0 + ''^ and consider negative values of t in the second equation, 
they become 

r- (2') 



u 



Dividing the first of these equations by the second, we get 
or, simply, 



r{^,-\-n).r{'^^-t) = R\ (28) 

The point ^-o we have termed a middle point, and B is the middle distance. 



(29) 
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Let us put 

4'^.= — C + G)^, 
4'. = — C + 6)„ . 

These points are also middle points. To show this let us put in (18) d = il/^-\-r. 
We get 

^(^, + T)= ^i-J^{ «. + io,,+ r), 

r r R^ 

-(^,+ r).-(^,-r) = ^, 

f (^, + t) . f (4^. - <r) = + -^, (30) 

'0 '0 '0 

The addition of 2a^, 2a^ or 2g)„ to any of these points brings us to another 
middle point. In the complete parallelogram of periods there are thus 16 middle 
points. 

Putting -r = i«;^ — ^ti, in (28), we get 

B,B^ = JR^, (31) 

which may be seen directly from equations (20). We find also, by putting 

X (32) 

Vy 
38 



Dividing, we find 
and similarly 
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§ 5. Reciprocal Orbits. 

If the initial conditions (equation (6)) be changed, that is, we consider 
another orbit which has the initial conditions 



d^ 



r, 







e = 0, r = ro, ~^ = — a 



(33) 



the constant of areas, h, being the same as before, and we determine r^ so that 

J^5 = iA l + a'-J^,=H (34) 

(6 and /? having the same quantitative values as before), and the differential 
equation of the orbit becomes 

''^- •-■ •-■' (35) 



?=*''-© +*Kf)'- 



de / ^r^y \ro 

This equation is the same as (8), except that b and /3 are interchanged. 

It is indeed the same as equation (9). The solutions for — and its reciprocal — 

are obtained at once from (18) by interchanging b and ^ (e^, e^ and e^ re- 
maining unchanged) and changing the signs of the solutions on account of the 
chans^e in sign of the derivative. We have then 

|= + ^l^^(5 + ^-i"-). 

From the second of (36) and the first of (18) we have, by virtue of the initial 
conditions. 

Therefore 
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Substituting this value in the first equation of (36) it becomes 



n 



>||-^(^ + o-4"J (38) 



Comparing this with the second equation of (18), we have, when i = 0, 



and consequently. 



or. 



We will call this solution the reciprocal solution. 

If in equation (37) we replace /"o hy ^o, using the formula 



r 

»-o 


u 


z:: 


. ^0 

r 


rrz 


= roro 


= 


■.B\ 


r 


f 


= 


1. 



it becomes 



(39) 



'''o — ^ -g- *'o> 

= +a!|^^(^ + ' + *".+ 2"a), (40) 

Comparing this with the first equation of (18), we find 

r (0 + G)^ + 2(0,) = r (0 + 6), + 2c),) = f (Q). (41) 

The reciprocal solution is thus included as a branch of the original solution. 

§ 6. Velocity in Orbit. 
The expression for the velocity in orbit is given by 

which is the same as equation (4). The value of Cj, as there determined by the 
initial conditions, is 
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Substituting this value of Cj , we have 

«^=i&(|j + ^). (42) 

At the middle point t^qj 

§ 7. Tangents and Asymptotes. 

Denoting by 4> the angle between the radius vector and the tangent to the 
curve, we have from the calculus 

tan^ = ^. 

dd 
Hence 

which reduces without difficulty to 

i6r§r^[-l,+ i:]8in2<^ = l, (44) 

or, as may also be written, 

r^ v^ sin^ ^ = ¥, 



Wy[J. + ^] 



sin^<^ = 1. 



Since r sin ^ = j) is the perpendicular from the origin to the tangent, we have 



h 

V = — . 

P 



The derivative of tan ^ with respect to r vanishes for ( — j = -5^; that is, 

at the middle points. 

The expression for the polar subtangent, 



hy a Force Varying Inversely as the Fifth Power of the Distance. 293 

r , Tn 

has, as — ^ « , the finite limiting value . ° , or, which is the same thing, 

— f-=^. This limiting value of the polar subtangent is the perpendicular 

distance from the origin to the asymptote. Consequently every infinite branch 
has an asymptote. 

§8. Determination of the Time. 

From the integral of areas, 

"^ dt "' 

we have, by substituting the value of r^ from (17), 

J [F(0 + c + ia.,)-eJd6=A.d<, (45) 

which can be integrated directly, since the o function is the negative derivative 
of the ^ function. 

Integrating and determining the constant so that Q vanishes with t, we get 

^(c + io,) -a^ + « + 4"J - ^-^^ = 0<- (46) 

Consider now another orbit (the symbols for which we will denote by the 
subscripts 1) such that 

A ro = hiTi. 
Then we will have also 

r = Yu 

R Ri 

Under these conditions the orbits will be similar, but will differ in size. The 
expression for the time will be given by 

? (ci + i CO,) - ^ (0, + c, + J 6.,) _ e, 0, = 1^ t,. (48) 
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Now if we take 6 = di and compare (46) and (48), we find 

which by virtue of (47) becomes 

£ _ r| _ ^ 

We have thus the following analogue of Kepler's Harmonic Law : 

Theorem : Corresponding arcs of similar orbits are described in times which 
are proportional to the cubes of the middle distances. * 

Kepler's Law does not depend upon the eccentricity or shape of the orbit, 
while the present law is restricted to similar orbits. 

Part II. 
§ 9. Classification of Real Orbits. 

It has already been observed that e^, e„ and e, are unaltered when b and (3 
are varied in such a manner as to keep their product b^ constant. If we put 
5/3 = y, then for fixed values of y different values of b correspond to different 
starting points on the orbit. Aside from questions of scale and orientation, there- 
fore, we find all the essentially distinct orbits as a one-parameter (y) set of curves. 
Restricting ourselves hereafter to realvalices of the variables, these orbits fall naturally 
into two classes : First, orbits having apses ; second, orbits not having apses. 

The condition for an apse is 

dS "' 
and from (8) we see that this condition is attained when — passes through a root of 

*(i) = "-(0 + *K^)* = ''- 

If, however, these roots are all complex, the condition for an apse can not be 

satisfied for any real value of — . Considered as a quadratic in ( — j , R f — j 

will have complex roots if its discriminant is negative; that is, if 

1 — 5/?<0. 

*Stader gives this result for circular orbits abont the origin as center. 
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From this condition we see that if h(i -^y is less than 1, the orbit will have an 
apse, and for values of y greater than 1 there will be no apse. Expressed in 
terms of y, 

e. = * + J ^T^, 

e.= -i, (49) 

e, = i — J V 1 — y. 

For values of y <^ 1 these expressions are all real. For values of y >• 1 e^i 
and e^ are conjugate imaginaries. 

When 6),, e^ and e are all real, it is customary to assign the subscripts in 
such a manner that 

For all values of y between — wand we see from (49) that 

e^ > e^ > e, , 

and consequently in this range of y 

;i = 1, |w = 2, V = 3. 

When y lies between and + 1> 

«A > e. > e^ , 
and in this case 

W= 1, V = 2, (i = Z. 

When y is greater than -\- 1, e^ and e, are complex, and we follow Weierstrass 

in taking 

;\, = 1, itf = 2, r = 3 (discriminant negative). 

We have therefore three distinct cases depending upon the value of y. 
There are also two limiting cases, y = and y = 1, the solution for which can 
be expressed by means of trigonometric and logarithmic functions. 

§ 10. Case I. — 00 < y < 0. 

We have already mentioned that in this case e^, e^ and e„ are real and that 
the subscripts have the order 

X = l, /M = 2, v = 3. 

With these values equations (17) become 

(Ty=itp(e+«+i"2)-^j, 
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This last equation is what we have called the reciprocal solution (see equations 
(33) et seq.) 

Since y = &/? is negative, either /3 is positive and h negative or the reverse. 
If |5 is positive and c = Ui — iog, then at = the particle is at an apse. This 
value of c, substituted in (50), determines the corresponding values of h and /? 
to be X Q 



h =63 — 63 



"With these values equations (50) become 

P {e + "1) — 63], 

(51) 






and consequently 



( 
(0 = 



63- 



[p(0 + «3) — ^], 



r 



\ 9{B)—e^ -y_ 0-3 ,.. 
\.^iB)-eA- c^'^' 

Since h = T3— 3 and is negative and If and t% are certainly positive, h must be 

negative. That is, the particle is moving under a repulsive force in the first 

h 
equation of (52). ^ = j^^ is positive. Therefore Te is positive and the particle 

is moving under an attractive force in the second equation. If the signs of b 

and /3 are interchanged, . the value of c at an apse is c = coi -^ ia^, and the 

solutions (52) are merely interchanged. 

The real half period % and the purely imaginary half period 03 are given 

by the formulae* 

E iK' 

"1 = / , "3 



V- > "'■6 / = 1 

gj — 63 V 61 — 63 

where ^ ^ 

K= r # K' = r # 

Jo Vl— ^sin^4)' Jo Vl— i'^sin^^' 

and ^ _ e^ — % _ V" 1 — y — 1 



Ci — 63 2 vT^ 



;fc'2 = 



_ 61 — eg _ V" 1 — y + 1 



Ci — 63 2^1 



* Schwarz, § 37, (3). 
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As B increases from to ui, p(0) decreases from + <» to Cj. Consequently 

r r 

— increases from 1 to -f <» and — decreases from 1 to 0. For the computation 

of intermediate values the 3 functions of Jacobi are convenient, and we find* 
where v = - — and r = ^ . 

2 6)x % 

The expression for the time (46) for — , as measured from the apse, becomes 

i0 + a"i)-r(^ + %) = Hi + n/T^) ^ t. 

By differentiating logarithmically the formula f 



a,(e) = e-'K*?L(i+^), 



where riy = ^ (o„), we obtain 



^ (0) = _ ,7, + ?^ (0 + ca.) = ^e + o.) - a^.). (54) 



Hl + Vl-y)^^=i0-SW- (55) 



Putting r = 1 in this expression, the above expression for the time may be 
written 

This expression for t vanishes for 6 = and becomes infinite for = Uj . 

The formula for the time for - is obtained in the same manner as the above 

differing only by the interchange of the subscripts 1 and 3. Making these 
changes, we get 

i(v/T^r^-i)^=P-^(0)-i0. (56) 

As 6 increases from to w,, t increases from to mJUX giJ — ^ at which 

A(Vl — y— 1) "' 

time the particle arrives at the origin. 



•Schwarz, §45, (11), (13). tSchwarz, §18, (1). 

39 
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§11. Case 11. 0<y< 1. 

In this case e;^ >• e, >■ e^ , and therefore 

% — l, v = 2, 11 = 3. 
The derivation of the formulae is the same as in the former case, the only difference 



throughout the work being the interchange of the subscripts 2 and 3. We find then 

where 

K iK' 



Vei — eg' Vei — eg' 

and 

— ^2 — ^— 1 — ■V^ 1— y 7V2— 2V"T 



]^ — f? 5 ^^ i -y 1 y ^2 _ 



(57) 



«! — «8 l + \^l— y' 1 + Vl— y* 

In terms of the S^ functions these solutions are 

f =^(e) = VT^(r|r). (68) 

'0 ^1 ~2 

The expressions for the time are 

(fort) 4(1-VT^)| =i«-iW- 

§12. Case III. !</< + <». 
For this range of values of y we will use the notation 



(59) 



^1= i + ii V y — 1, 
/> — 1 

H Ti } 

63 = 1 — JiVy — 1; 
so that ^ = 1, |u = 2, v=3 (discriminant negative). With this ordering of the 
subscripts the half periods wj and o^ are conjugate imaginaries, while (aj^Oi+Og 
is real. For real values of the argument the p function varies from -f oo to ^2 
as the* argument varies from to Og. 
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The solution (17) then is 

(fy=l-[p(e+«+*"2)-«2], 

The apses in this case are all imaginary, but the middle point i^o *°<i the 
zero-point B^ are real. If = at the middle point, then c = 0; and since 
j? (J ©2) = p ( — J02), we see by comparing the two equations (60) that 
& = |S = V y, and therefore 



(61) 



therefore 






kJ y a 






(62) 



V y <r 
and comparing these two equations it is evident that 

(«) = 



or 

r(e).r(-0) = ri. 

The formula (43) for the angle between the radius vector and the tangent to the 
curve becomes in this case 

Since ( — } and its reciprocal enter this expression symmetrically,^ has the 

same value for + 6 as for — %. It is also readily verified that ^ is a maximum 
for = 0. 



(65) 
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If we put 

^ — e3 = pe«*, p>0, 

eg — Ci = pe-**, 0<'4/<7r, 

we find the following value for og • * 

Vp 
where kl = sin* i-^*. 

Solving (65) for p and 4', we find 

p = J V y, ^ 

sin='J,^=l±^. 
2 V y 

The formula for the time reduces to 

^^^'r^i="^^-^{^ + i<^) + ^{i^z)^ (66) 

»o 

§13. Limiting Gases. 
a) Limiting Case y = 0. 

Since y = J /3 = 0, either 6 or /3 is zero ; and since the solution for either 
is the reciprocal of the solution for the other, it is immaterial which we choose. 
We will take (3 = 0, and then it follows from their definition (7) that 

The differential equation (8) becomes 

The solution of this equation is 

^ = VT^ cos (6 + c). (68) 

'0 

If = at an apse, we find c := and consequently J 5 = i, and the solution 
becomes 

— = cos 6, (69) 

* See Appell et Lacour, "■ Thiorie des Fonctiont Elliptiqves," §138. 
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which is a circle passing through the initial point and the origin and having this 
line as a diameter. 

The reciprocal solution is 

^ = sec e, (70) 

which is a straight line and implies that the force acting is zero. 

b) Limiting Oase y = !• 
Since b^ = 1, b= -n and the differential equation (8) becomes 






± 



^[0-*]- <"' 



de V 

The solution of this equation is 

where, to satisfy the initial conditions, we must take 

c=logl-^j. 
1 + V J 

This solution, together with its reciprocal, represents two spirals asymptotic 
to the circle r = y^ V b, one lying within the circle, the other without the circle. 
For the special values 5 := /? = 1 this solution degenerates into the asymptotic 
circle itself; that is, 

r = rp. 

/r \^ dd h 
From the relation ( — ) -^j = -j and equation (20) we derive the following 

expression for the time in terms of 6 : 

h _ 1 4_ />V2» + c 

^,(,_« = »+V2l±i,^. (73) 
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§ 14. Eesumi. 
Tablb of Cases. 



r 


«i 


e. 


e. 


U) 


*' 






-i 
-i 
i-i Vi-r 
-i 
-i 




K 

"'•"vr-r 

:r/2 

2ir 
' 1+Vi-r 

+ 00 






V 1-r-i 
2 V i-r 




r<0 


t+iVi-r 


i-J V 1-r 

-i 

-i 

i 




1- Vi-r 
1+ Vi-r 

1 

1 + V7 
2V7 


o<r<i 


i+ J Vi-r 
* 


r>i 


i + it Vr-i 



Table of Ebsults. 



\l 



ifi 



r<o 






(«, — ei) 









= 



cos e 
sec B 




1 



o<r<i 



= 1 






(«2 - «») 
(«i - «») 



?'>1 



Vr 



V r " 



i 



}V" 






^[J^ + isinae] 



rj 



tan (7 



A(«i-e»)L -^i J 



i 






iV" 



2r5 



^[i<>-C(J<», + ^) + C(Ja,,)] 



t>(0)=+oo, !?(««,) = e,, l?K) = e,, ?(<«,) = «,. 
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The expressions in this table assume the initial position at an apse or at a 
middle point. 

§15. Orbit Considered as Fwnction of Parameter y. 

It is interesting to trace the changes in the form of the orbit as the 
parameter y runs through its range of values from — oo to + « . For 
Y ■=■ — 00 , 0)1 = and the orbit is a straight line through the initial point 
and the origin. 

For finite negative values of y there are two orbits, one lying wholly within 

the circle, — = 1, the other lying wholly without the same circle (Fig. 1), and 




y = -i- 



Fig. 1. 



-f =^W. 






SO related that the outer orbit is the transform of the inner by reciprocal radii. 
The inner orbit consists of a series of loops passing through the origin and 
tangent to the circle, repeating themselves at intervals of 2 Oj . The curve is 
closed only when o^ is commensurable with n. 

As Y increases towards 0, Oi approaches J n. The loops of the inner orbit 
broaden out and approach the circle having the initial point and the origin as a 
diameter. The outer orbit expands very rapidly and approaches the straight 
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line tangent to the circle at the initial point. These limits are attained for 
ly =z (Fig. 2). These outer orbits are denoted by (a) in the table. It will be 
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observed that the corresponding value of ib is e^ — e^, which is negative. But 

by definition J 6 = — ^g — g , which for real orbits can be negative only if Je is 
2 fi rQ 

negative. This means that the outer orbits are described under a repulsive force 

instead of an attractive force. 

As y passes through the value of b for the outer orbits also passes 
through and becomes positive; the force is therefore changed to attraction. 
It will be observed that for these orbits the expression for the radius of curvature 
also changes sign. The orbits therefore curve in the opposite direction. As y 
increases from to 1, the curvature increases, oj increases towards + oo , and as 
y =1 the orbit approaches the hyperbolic spiral (equation (72)) asymptotic to the 
circle. The loops of the inner orbit continue to widen out (Fig. 3), and as y = 1 
the inner orbit also approaches the hyperbolic spiral asymptotic to the circle on 
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the inside. These limiting spirals are actually attained for the value j/ = 1 if 
the initial position is not taken at an apse. Otherwise the limiting orbit is the 




Vui. 4. 

/ = 1- 



asymptotic circle itself (Fig. 4). In order to avoid confusion in the figure only 
that branch of the outer curve has been drawn corresponding to positive values 




201 

''= 200 • 



FifS. 5. 



(l6)„ + f)). 



of d, and of the inner orbit the branch corresponding to negative values. As y 
passes through the value 1 these branches which are asymptotic to the circle 
40 
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unite, and for values of y >■ 1 the orbit crosses the circle and the apse ceases to 
exist. The point of crossing becomes a middle point, and the reciprocal relation 
between the inner portion and rhe outer with respect to this point continues as 




7 = 16 : 



i'lG. 6. 
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vy 



■ (J u, + w. 



before. For values of y but little greater than 1, the period o^ is very great, 
but decreases as y increases. The curve which now passes through the origin 
and infinity straightens out as y increases (Fig. 6). As j^ == -(- tao , ojg decreases 
and approaches the value 0, and the orbit again approaches the straight line. 

The University of Chicago, June 27, 1907. 



